Thickness dependence of surface diffusion in epitaxial LaAlOs on SrTiO3(001) 
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The LaAlO3/SrTiO3{001) thin fihn materials system was studied using in situ, simultaneous x-ray 
diffuse scattering and specular reflectivity during pulsed laser deposition. Using this method, we 
are able to measure the time dependence of the characteristic surface length scale and the charac- 
teristic time for both in-plane and downhill diffusion. These data allow for the determination of the 
activation energy for various diffusion processes as a function of LaAlOs thickness. Additionally, we 
show that the downhill diffusion rate of the first monolayer is distinctly different than subsequent 
layers. These results are directly compared to previous experimental observations seen during the 
deposition of homoepitaxial SrTiO3{001). 



Complex oxides exhibit a vast range of materials 
properties encompasing electrical insulators, high-Tj, su- 
perconductors, semiconductors, dielectrics, ferromagnet- 
ics, and multiferroics. Of particular note, these ma- 
terials admit the deterministic synthesis of essentially 
arbitrary crystalline structures, which is leading to 
new discoveries. For example, for certain thicknesses, 
the LaAlO3/SrTiO3(001) system exhibits a conducting, 
quasi-two-dimensional electron gas at the interface be- 
tween two wide-bandgap insulators [IHl]. The forma- 
tion of the conducting layer is sensitive to the number of 
LaAlOa layers [5], which leads to the possibility of new 
device architectures [6]. Pulsed laser deposition (PLD), 
long known as the tool of choice in a research context for 
the growth of complex-oxide materials, has now demon- 
strated the ability to grow interfaces with a sharpness 
previously thought attainable only with molecular beam 
epitaxy. With the keen interest in growing such het- 
erostructures and artificial materials via PLD, it is crit- 
ically important to understand the microscopic growth 
processes setting the fundamental limits to obtaining 
atomically perfect growth. Here, time-resolved, simulta- 
neous diffuse and specular x-ray scattering measurements 
are used to determine the rate of several surface diffusion 
processes during the PLD of LaAlOs on SrTiO3(001). 
Perhaps not surprisingly, the phenomenology we observe 
in this heteroepitaxial system is much richer than that 
observed in homoepitaxy, with film thickness dependent 
kinetic coefficients. We find that the barrier for in-plane 
diffusion of LaAlOa on SrTiOa is substantially larger 
than the rate-limiting barrier for downhill transport from 
an LaAlOa island to the SrTiOa substrate. And, the en- 
ergy barrier for diffusion is found to be a function of 
LaAlOa film thickness, for films from one to six unit cells 
thick. 

These studies were performed in the PLD/x-ray 
diffraction chamber in the G3 experimental station of the 
Cornell High Energy Synchrotron Source (CHESS). The 
supplemental documents of Ref. [7] contain the experi- 



mental details. A one-dimensional diode array detector 
(pixel size 0.125 x 5 mm) was used to acquire the dif- 
fuse scattering data at the quarter-Bragg position on the 
crystal truncation rod. A KrF excimer laser (248 nm), 
focused down to a 7.4 mm^ spot, was used to ablate the 
single crystal LaAlOs target. The laser fluence was 1.6 J 
cm^2 for all films, producing a deposition flux of approx- 
imately twelve pulses-per-monolayer. Here, a monolayer 
(ML) refers to a one unit cell thick layer. The laser was 
fired at a repetition rate of 0.19 Hz. The substrate was 
held at controlled temperatue in a background pressure 
of 7.5 X 10~^ Torr of O2. All films were grown in a layer- 
by-layer growth mode to a thickness of six MLs, which is 
well below the critical thickness for strain relaxation [8]. 

In layer-by-layer growth, islands nucleate, grow, and 
then coalesce. At high coverage, near layer completion, 
the surface is best described by a series of pits that fill in 
as more material is deposited. During deposition, a char- 
acteristic length scale Ligi, is present on the surface. This 
length scale is a measure of the average separation be- 
tween islands (or pits) ^9J. At low coverage, x-ray diffuse 
scattering directly measures both Lisi and the character- 
istic time for in-plane (intra-layer) mass transfer, Tisi [7] . 
By monitoring the surface roughness using the specular 
reflectivity, the characteristic diffusion time for downhill 
(inter-layer) transport, Tspec, may be determined [71 1101 - 
[T^ . Identifying the diffusion length as, Ld = Lisi/2, the 
intra- and inter-layer diffusion rates may be calculated 
by invoking the Einstein relation: D ~ L^/^t jl4j . 

The x-ray scattering may be broken into specular 
/specand diffuse Idiff components and each component 
can be accurately described by the functional form 

I{q\\) = h/[l + e{%-q,?f'\ (1) 

where q\\ is the in-plane scattering vector, and 50 is the 
peak position [7]. For Idiff, Qo is directly related to 
the distance between islands or pits by the equation: 
qq « 2iT /Lisi [15]. For /spec, <Zo = 0. During homoepitax- 
ial layer-by-layer deposition, Ispec oscillates due to the 
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FIG. 1; (a) X-ray scattering data for the PLD of LaAlOa 
on SrTiO3(001) at 781°C. The intensity measured at / 
is the diffuse x-ray scattering. The diffuse scattering is 
directly related to the LaAlOs island distribution, (b) The 
corresponding fit results for Ispec (black) and lisi (green). 
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FIG. 2: Diffuse scattering data and the corresponding fit 
results for the first ML of LaAlOs at 1000° C and 858° C. The 
vertical lines represent the times when the laser was fired, (a), 
(c) and (e) correspond to the 1000° C growth and (b), (d) and 
(f) show the 858° C results. The arrows annotate the data 
points used to determine nai and Tapec- 



continuous roughening and smoothening of the surface. 
Each local intensity maximum corresponds to the com- 
pletion of wl ML [71[T0l[12]. In heteroepitaxy, interfer- 
ence between beams reflected by the interface and the 
surface create an additional modulation. These inter- 
ference effects are known as "Kiessig" fringes, and their 
periodicity is dependent on the out-of-plane scattering 
vector, |16j . For our scattering geometry, the enve- 
lope of Ispec has a periodicity of ~4 MLs, and the rough- 
ness oscillations within this envelope have maxima cor- 
responding to the completion of wl ML. 

To interpret the data, we write the total scattering in 
a Qz plane of reciprocal space as 

Itot Ispec ^" lisl ^" Ism: (^) 

where lisi is the observed diffuse scattering [21], and Ism 
accounts for scattering by features smaller than our de- 
tection limit. (Since our maximum measurable value 
of (7|| is Qmax ~ 0.2 A~^, wc are insensitive to surface 
fluctuations with wavelengths smaller than 3lA.) Con- 
sequently, if the coverage is low and I spec is constant, a 
time dependent lisi results from mass transfer between 
very small islands to islands measurable by our experi- 
ment. Therefore, the time dependence of lisi provides 
a direct measure of the time scale for intra-layer mass 
transport T]. 

Figure flla) shows a false color image of the scattered 



x-ray intensity, as a function of time and q\ | , during the 
deposition of six MLs of LaAlOs on SrTiOs. Each 260- 
msec time slice, represents intensity vs. (see Ref.[7] 
for details). The intensity seen at g|| 7^ is a result of 
the Henzler ring of diffuse scattering |17j and, as men- 
tioned above, the size and shape of this ring is deter- 
mined by the LaAlOa island distribution on the surface 
[9]. As observed during PLD of homoepitaxial SrTiOa 
[7], Qo decreases as each additional layer of the LaAlOa 
fllm is deposited, representing a systematic decrease in 
the island density with increasing layer number. To de- 
compose Ispec and lisi , each time slice is flt to our model; 
the results are shown in Fig. [ijb). In this flgure, lisi os- 
cillates with the period of 1 ML, and is out of phase with 
Ispec- As discussed above, both the roughness and Kies- 
sig components to Ispec are visible in Fig. [IJb). The 
periodicity seen in Ispec that is n out of phase with Iisi is 
the roughness component, and it results from the scatter- 
ing of x-rays by the LaAlOs islands ^18^ . The minimum 
in the envelope of Ispec at around 200 seconds is a signa- 
ture of the Kiessig component, therefore, the Bragg 
position sets a 4:1 ratio between the diffuse and Kiessig 
periods. 

To examine the influence of temperature on the 
LaAlOa deposition, we show false color images of the 
scattering data obtained during the deposition of the flrst 
ML at 1000° C and 858° C in Figs, ^a) and^h), respec- 
tively. Immediately noticeable in the data is the increase 
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in qo as the surface temperature is decreased. This im- 
phes that the island density increaseds with decreasing 
temperature, consistent with classical nucleation theory. 
These data were fit to our model, and the results are 
shown in Figs. [2|^c)-[2f^f). As seen in Fig. [2f^c), qo mono- 
tonically decreases both inter- and intra- pulse at 1000°C. 
This behavior is attributed to coarsening of the LaAlOa 
island distribution on the surface [H IH [19]. Fig. ^d) 
shows qo for the 858° C deposition. At, and below, this 
temperature, the low diffuse intensity precluded an ac- 
curate fit before the second laser pulse. An interesting 
feature of Fig. ^d) is the slight increase in qo just af- 
ter the third laser pulse. The best fit estimates for the 
peak positions just before and after the third laser pulses 
are qo = 0.078 ± O.OOlA^i and qo = 0.081 ± O.OOlA"!, 
respectively. This increase in qo, as additional material 
arrives indicates an increase in island density. Therefore, 
island nucleation may be occurring up to this coverage 
at 858° C. 

Figures [2]je) andpjjf ) show Igpec and lisi for the data 
presented in Figs. [2pi) and[2]jb), respectively. Following 
the third laser pulse, a constant Ispec accompanies the 
slow rise in lisi. As discussed previously, this temporal 
dependence of lisi corresponds to the timescale for in- 
plane mass transfer. Additionally, the time dependence 
of Ispec at high coverage is a direct measurement of the 
timescale for downhill diffusion. Fitting both lisi and 
Ispec to a simple exponential function allows for the de- 
termination of these characteristic diffusion times. 

Conceptually, the deposition of the first layer in het- 
eroepitaxy is distinctly different than deposition of subse- 
quent layers. For example, the diffusion rate for LaAlOa 
species on the SrTiOs substrate need not be the same 
as diffusing LaAlOs species on the LaAlOa film. Ad- 
ditionally, a mobile species diffusing on top of the first 
LaAlOa ML might interact with the SrTiOa substrate 
during downhill diffusion. Therefore, one might expect a 
film thickness dependence of the diffusion rate until the 
LaAlOa film reaches a critical thickness for substrate in- 
teraction effects. To examine this possibility, the inset 
to Fig. [s] shows a plot of vs. Tspec for each of the 
6 ML deposited at two substrate temperatures: 606°C 
and 704°C. We note that after the first ML, and Tgpec 
show a linear relationship, where the slope of the line is 
related to the diffusion rate by: — ADTgpec- The fact 
that the first data point is inconsistent with the remain- 
ing points suggests that the interlayer transport for the 
first layer is distinctly different than that of subsequent 
layers. 

To further examine the diffusion behavior in the 
LaAlOa films, an Arrhenius plot is shown in Fig. [3] for 
three data sets: the intra-layer diffusion rate for the first 
ML, the inter-layer diffusion rate for the first ML, and 
the average of the inter-layer diffusion rates in the sec- 
ond through sixth MLs. The intra-layer diffusion bar- 
rier for LaAlOa species on SrTiOa was determined to be 



Eai = 0.88 ± 0.33 eV. The noise in lisi precluded the 
determination of Tisi for subsequent MLs, therefore, the 
intra-layer diffusion rate for LaAlOa species on LaAlOa 
could not be measured. Fitting the inter-layer diffusion 
data to an Arrhenius model yielded activation energies 
oiEa2 = 0.40 ±0.06 eV and = 0.19 ±0.10 eV for the 
first ML and second-sixth ML, respectively. These data 
illustrate the layer dependence of the activation barrier, 
showing that the largest energy barrier corresponds to 
the diffusion of LaAlOs species on the SrTiOa substrate. 

The various diffusion processes measured from our 
data are presented schematically in Fig. |4j The diffu- 
sion barrier, Eai, corresponds to in-plane mass transfer 
of LaAlOa species on the SrTiOa substrate. The activa- 
tion barrier, Ea2, corresponds to the inter-layer diffusion 
from the top of the first LaAlOs layer down to the sub- 
strate. i?a3 corresponds to the downhill diffusion of the 
second-sixth MLs to the underlying LaAlOs layer. Our 
measurements do not allow us to determine if the inter- 
layer activation barriers are the Ehrlich-Schwoebel (ES) 
barrier, the surface diffusion barrier, or the sum of the 
two. 

Our results explicitly demonstrate the contrast be- 
tween heteroepitaxial and homoepitaxial diffusion pro- 
cesses. Specifically, previous work measuring diffusion 
rates for PLD of homoepitaxial SrTiOa showed that the 
energy barriers for inter- and intra-layer diffusion were 
the same: Ea = 0.97 ± 0.07 eV [7]. It was also shown 
that the inter-layer diffusion rate is the same for all MLs 
during SrTiOa homoepitaxial deposition (see Fig. 3 of 
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FIG. 3: Arrhenius plot of the diffusion rate for LaAlOs 
species during three different processes: in-plane diffusion on 
the SrTiOa substrate (blue, triangles), downhill diffusion from 
the first LaAlOs ML to the substrate (red, points), and down- 
hill diffusion from a LaAlOs ML to the LaAlOs film (black, 
circles). The inset shows the dependence of Ld on r^pec for 
each ML at 704° C (squares) and 606° C (circles). The arrows 
point to the data points for the first MLs. 
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FIG. 4: Schematic description of the activation barriers mea- 
sured in Fig. [3] The oscillations show the various activation 
energies across the surface, with the large increase at the step 
edge representing the ES barrier. Eai is the in-plane diffusion 
barrier for LaAlOs species on the substrate. Ea2 is the en- 
ergy barrier for LaAlOa species diffusing down to the SrTiOa 
substrate. Ea3 represents the downhill activation barrier for 
layers 2-6. The dotted outline illustrates which set of diffusion 
processes determine Ea2 and Ea3- 



was used to explain the persistent layer-by-layer growth 
mode. Additionally, we have shown that all diffusion 
barriers for LaAlOa on SrTiOa are small compared to 
homoepitaxial SrTiOa. 
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ref. [7]). Here, the measured intra- layer diffusion barrier 
of the first LaAlOa ML is found to be larger than the 
downhill activation barrier. Additionally, the inter-layer 
diffusion of the first ML of LaAlOa on SrTiOs was found 
to differ from subsequent layers, illustrating a distinct 
difference between homoepitaxy and heteroepitaxy. 

Our data may be used to explain the persistent layer- 
by-layer growth mode for the LaAlOs/SrTiOa deposi- 
tions. To do this, it is useful to consider the Volmer- 
Weber (3D) growth mode. In Volmer- Weber growth, 
the lower diffusion barrier of the film/substrate inter- 
face, when compared to the film/film diffusion barrier, 
results in the sticking of adatoms to the top of nu- 
cleated islands. This behavior has been seen in Fes Si 
growth on GaAs(OOl), where the growth mode was di- 
rectly correlated to the surface diffusion barriers pO] . For 
LaAlOa/SrTiOs, Eai > Ea2, providing the proper con- 
dition for inter-layer transport, a prerequisite for layer- 
by-layer growth. Therefore, diffusing species will tend 
to absorb to the substrate rather than the film. Addi- 
tionally, this physical description may be applied to the 
deposition of the second ML of LaAlOs, since Ea2 > Eaz- 

In conclusion, we have shown that simultaneous, time- 
resolved surface diffuse and specular scattering can be 
used to measure the surface diffusivity during the pulsed 
laser deposition of LaAlOa on SrTiO3(001). We have 
used this data to measure energy barriers for three dif- 
ferent surface diffusion processes: in- plane diffusion of 
LaAlOa on SrTiOa, downhill diffusion of LaAlOa to the 
SrTiOs substrate, and downhill diffusion of LaAlOa to 
the LaA103 film. The activation barriers were found to 
decrease, respectively, for each of these processes, which 
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